Boundary instability of a two-dimensional electron fluid 
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It was shown previously that the current-carrying state of a Field Effect Transistor with asymmet- 
ric source and drain boundary conditions may become unstable against spontaneous generation of 
plasma waves [l[ . By extending the analysis to the two-dimensional case we find that the dominant 
instability modes correspond to waves propagating in the direction perpendicular to the current and 
localized near the boundaries. This new type of instability should result in plasma turbulency with 
a broad frequency spectrum. More generally, it is shown that a similar instability might exist, when 
a strong enough current goes through a single boundary between the gated and ungated regions. 



It was shown previously 1[ that the current-carrying 
state of a Field Effect Transistor may become unstable 
against spontaneous generation of plasma waves in the 
transistor channel, provided there is an asymmetry in the 
boundary conditions at the source and at the drain. An 
extreme case of such asymmetry is the ac short-curcuit 
condition at the source and the ac open curcuit at the 
drain. For submicron gate lengths the frequencies of the 
plasma oscillations belong to the terahertz range, thus 
the FET can, in principle, serve as a generator of tera- 
hertz radiation. The nonlinear properties of the electron 
fluid in the transistor channel can be also used for detec- 
tion and frequency mixing in the terahertz domain 0]. 

Experimentally, both terahertz emission [H, 0, Q and 
detection Q in nanometric transistors were demon- 
strated. Fig. 1 presents experimental data Q for a 
GaAlN/GaN HEMT at 4 K clearly showing the emission 
threshold at a certain source-drain voltage (or current) 
and a typical broad emission spectrum in the terahertz 
domain. Contrary to the prediction of Ref . [1] , the spec- 
trum depends neither on the gate length, nor on the gate 
voltage. Similar results for terahertz emission were ob- 
tained at room temperature [5[. 

It is not firmly established that the observed emission 
is indeed related to the instability predicted in |l| (see 
5]). However, one cannot directly compare the theory 
with the experiments because the experimental geometry 
is very different from the one-dimensional model adopted 
in [l[ . In the standard experimental situation, the width 
of the gate W is much larger than the gate length L, 
typically W/L ~ 100, see Fig. 2, left. Under such condi- 
tions the one-dimensional model, where the plasma den- 
sity and velocity depend on the coordinate x only, is not 
appropriate, since obviously oblique plasma waves with a 
non-zero component of the wave vector in the y direction 
can propagate. In such a geometry, the gated region is 
not a resonator, but rather a waveguide with a continu- 
ous spectrum of plasma waves, see Fig. 2, right. 

The purpose of this work is to extend the analysis of 
stability of the steady-state flow [l[ to the more realistic 
geometry of Fig. 2. Since W » L, we will consider the 
limit of a strip, which is infinite in the y direction. It will 
be demonstrated that in such a geometry a new mode 
of instability dominates, which is localized near the gate 
boundaries. Moreover, a similar instability should exist 
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FIG. 1: Experimental results for THz emission from a Al- 
GaN/GaN HEMT at 4.2 K 0. a- The drain current (right 
scale) and the emission intensity (left scale), as functions of 
the source-drain voltage, U 3 d- Note the pronounced threshold 
for emission, b - The emission spectrum at U s d = 3 V. 



in the limit L — > oo, i.e. near a single boundary of a 
current-carrying two-dimensional plasma. 

Within the hydrodynamic approach the electrons in 
a gated 2D channel can be described by the following 
equations [l|: 



dV 

~dt 



(VV)V 
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FIG. 2: Left panel: the geometry of the gate. The width W 
is much greater than the length L. Right panel: The plasma 
wave spectrum in a strip, z = 2Lu/s is the dimensionless 
frequency, p = 2qL / s is the dimensionless wave vector in the 
y direction, s is the plasma wave velocity. 



dU 

at 



V(E/V) = 0, 



where V(r, t) is the electron hydrodynamic velocity, 
U(r, t) is the gate-to-channel voltage swing, r is the vec- 
tor in the 2D plane, e and m are the electron charge and 
the effective mass respectively. Equation (1) is the Euler 
equation, and Eq.(2) is, in fact, the continuity equation 
since the electron density in the channel, n, is related to 
the voltage swing, U, by the relation 



en = CU, 



(3) 



where C is the gate to channel capacitance per unit area. 
This equation holds if the scale of the variation of the 
potential in the channel is large compared to the gate-to- 
channel separation d (the graduate channel approxima- 
tion). 

Collision processes give an additional term — V/t in 
the right-hand side of Eq.(l), where r is the momen- 
tum relaxation time. In the following, this term will be 
neglected, however it should be understood that the in- 
stabilities studied below will practically exist only if the 
instability increment is greater than 1 /r, a condition that 
determines the instability threshold for the drift velocity, 
similar to the situation in the one-dimensinal model 

We chose the x axis in the direction from source (x = 0) 
to drain (a; = L) and, following Ref. 1, we impose the 
asymmetric boundary conditions of a fixed voltage at the 
source and a fixed current at the drain: U = Uo at x = 
and j x = jo at x = L, where j x is the x component of the 
current density. Because of Eq. (3), the latter condition 
can be rewritten as (UV x ) x =l — UqVq, where vq = jo/en 
is the electron drift velocity. 

It was pointed out in [l| that Eqs. (1) and (2) are 
identical to the equations describing the so called "shal- 
low water" in conventional hydrodynamics Q, plasma 



waves in the channel being analogous to water waves in 
the case when the wavelength is much larger than the wa- 
ter depth. Furthermore, it was shown that the current- 
carrying steady state described by the stationary solu- 
tion of Eqs. (1,2) with the above boundary conditions, 
U = Uo, V x = Vq, is unstable against spontaneous gen- 
eration of plasma waves with a growth increment given 
by 



7 = 



2sL 



■In 



vo i 



S - Vq' 



(4) 



where s — (eUo/m) 1 / 2 is the plasma wave velocity. 

This result followed from a one-dimensional analysis, 
e.g. small perturbations of the steady state were assumed 
to be independent of the coordinate y in the direction 
perpendicular to the current. As we shall see, the ex- 
tension of the analysis to y-dependent perturbations not 
only gives corrections to Eq. (4) but, somewhat unex- 
pectedly, gives a new mode of instability which always 
dominates. 

We study the time dependence of small perturba- 
tions of the steady state. Accordingly we put U — 



(2) Uq + (m/e)u,V x = v + v x ,V y 



v y and we linearize 
Eqs. (1,2) with respect to the small quantities u, v x , v y . 
The boundary conditions become: 



^=0 = 0, {v u + s 2 v x ) x= 



0, 



(5) 



(zero ac voltage at the source and zero ac current at 
the drain). We look for the solutions of the linearized 
equations with u, v x ,v y ~ exp(— iwt + ikx + iqy), where 
k and q are the components of the wave vector in the x 
and y directions respectively. This procedure gives 



(u) — kvo)v x — ku = 0, 



(oj — kvo)v y — qu = 0, 



(6) 



(7) 



[uo — kvo)u — s 2 (kv x + qv y ) = 0. (8) 

The dispersion relation for the plasma waves follows: 

(u-kv ) 2 = s 2 {k 2 + q 2 ), (9) 

the term kvo taking into account the Doppler shift due 
to the motion of the electron fluid. For given u> and q we 
find two values for the ^-component of the wave vector, 
corresponding to oblique waves propagating downstream 
and upstream: 



h,2 = 



-ljvq ± s^Juj 2 — (s 2 — v 2 )q 2 



(10) 



For q — this reduces to h\ = lu/(s + vq), ki — —lo/ (s — 
vq). The general solution for u and v x can be found using 
Eq. (6): 



Aexp(ik\x) + B exp(ilt2x) 



(11) 
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Aexp(ikix)-\ : Bexp(ik 2 x) 7 (12) 



u) — kivo 



u - k 2 v 



where A and B are constants. 

The boundary conditions, Eq. (5), together with Eq. 
(10) give the relation 



exp(i(fc x - k 2 )L) 



lu — kiVo 
u> — k 2 v ' 



which can be rewritten in the form: 



<■ m 2^ f3y/z 2 - p 2 - z 

cxp^Vz 2 - P 2 ) = - , = , 

P\J z z — p z + z 



(13) 



(14) 



where the dimensionless variables for frequency, wave 
vector, and drift velocity are introduced: 
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(15) 



Equations (14), (15) define the complex frequency lu — 
to + 17 as a function of the drift velocity Vq and the 
wave vector q. For q = one obtains the previous one- 
dimensional result [l| with lu' = irl(s 2 — v 2 ,)/ (2sL), where 
I is an odd integer, and the increment lu" = 7 given by 
Eq. (4). 

In the general case Eq. (14) can be solved only nu- 
merically. However, an analytical solution can be ob- 
tained for drift velocities small compared to the plasma 
wave velocity ((3 << 1). For (3 = the solution of 
Eq. (14) is z = (I 2 +P 2 ) 1 / 2 , or in dimensional units 

lu = s((nl/L) 2 + q 2 ) , which represents the spectrum 
of plasma waves in an infinite strip with the assumed 
boundary conditions at x = and x — L (Fig. 2). The 
linear in (3 correction to this value is purely imaginary: 
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('0 
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L 1 + {qL/(«l)f 



(16) 



Thus, as q increases and becomes comparable to or 
larger than the quantized value of k = irl/L for the l-th 
mode, the instability increment decreases from its value 
vq/L given by Eq. (4) for vq << s. The correction to 
the real part of lu is of second order in (3. 

However, in addition to this predictable result, another 
solution of Eq. (14) exists, for which z (or lu) is purely 
imaginary. For vq << s this solution can be found ana- 
lytically by assuming that \z\ « p and (z 2 —p 2 ) 1 ^ 2 w ip. 
This gives z = i/3ptanh(p/2) or, in dimensional units 



= and 



7 = qv<j tanh(gL). 



(17) 



For large qL this gives 7 = \q\vo, thus in contrast to 
the result given by Eq. (16), the growth increment of 
this new mode increases at large q, so that this mode of 
instability is the dominant one. The numerical solutions 
of Eq. (14) for (3 = 0.5 are presented in the Fig. 3, 
together with the approximate result for (3 « 1 given 
by Eq (17). 




FIG. 3: The instability increment z" as a function of the 
transverse component of the wave vector p in dimensionless 
units for j3 = vq/s = 0.5 (Eqs. (14, 15), numerical calcu- 
lation). 1, 2, 3 - for normal modes with I = 1, 3, and 5, 
respectively, 4 ~ f° r the new mode of instability, 5 - approx- 
imation given by Eq. (17). 



It can be seen that for qL >> 1 the new mode is 
localized near the boundaries at x = and x — L on a 
distance ~ l/q. For example, in the case (3 << l,qL >> 
1 we have fci ~ — k 2 ~ iq (see Eqs. (9,10)). Thus the 
instability mode is formed by waves which are evanescent 
in the x direction. 

Since for large \q\L the growth increment for the new 
mode does not depend on L, and since in this case the 
mode is localized near the boundaries, it seems plausible 
that a similar instability of the steady-state flow should 
exist for a single boundary of an infinite (both in the y 
and the x directions) two dimensional current-carrying 
plasma. We now show that this is indeed the case. 

Let a steady current with the drift velocity vq flow 
across the boundary (x — 0) of a semi-infinite sample 
situated at x > 0. The general boundary condition at 
x = is defined by the impedance C relating the ac volt- 
age and the ac current (compare with Eq. (5)): 



C(v u + s 2 v x ). 



(18) 



The boundary condition at x = 00 corresponds to the 
vanishing of the small perturbations, u = v x = v y = 0. 

The impedance C will be considered as purely imagi- 
nary: C = iX/s, where A is the dimensionless parameter 
proportional to the effective capacitance. (The existence 
of a finite resistance, described by the real part of £, 
will obviously introduce damping of the initial perturba- 
tions and, if it is large enough, any instability will be 
supressed). 

To insure the boundary condition at x — 00, we now 
keep only one exponent in Eqs. (11,12), with the wave 
vector fc, whose imaginary part is positive. These equa- 
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FIG. 4: The coefficient G in Eq. (20) for p = v /s = 0.5 
as a function of the dimensionless parameter A denning the 
boundary condition at x = (the boundary impedance is 
presented as £ = iX/s). For large A, G ~ 1/A 3 . 

tions, together with Eq. (18), give: 

iX 

lu — kvo = ask, a = — . (19) 

1 — iXp 

Inserting this relation in Eq. (9), we find the value of 
the wave vector k = i\q\(l — a 2 ) -1 / 2 , where the sign of 
the square root should be chosen so that its real part be 
positive. 

Finally, from Eq. (19) one finds ui. Its imaginary part 
7 defines the instability increment: 

7 = GkK G=^{0=). (20) 

Note, that the values of q in Eqs. (17,20) are limited 
by the condition q < 1 j d, where d is the gate-to-channel 
separation. For larger q the graduate channel approxi- 
mation, used in deriving Eqs. (1,2) breaks down. 

The dimensionless coefficient G depends on the value 
of A, defining the boundary impedance, and on the flow 
velocity Vq, see Eq. (19). For A << 1, we have G = 1 
and Eq. (20) coincides with Eq. (17) for large \q\L. With 
increasing A the coefficient G decreases (see Fig. 4), re- 
ducing the instability increment, which however remains 
always positive. 

Thus, if the condition \q\vo > 1/r is satisfied, the 
current-carrying steady state is unstable against small 



perturbations, and the region of instability is localized 
near the boundary. This is similar to what one observes 
in a river, when the water flows with sufficient velocity 
across an abrupt step in the waterbed: waves with wave 
vectors perpendicular to the flow are excited, while the 
wave vectors in the direction of the flow are purely imagi- 
nary, which accounts for the localization of the turbulent 
region near the step. 

Certainly, the linear theory cannot predict the out- 
come of this instability. However, since the spectrum of 
plasma waves is continuous, it seems likely that the in- 
stability will result in a turbulent motion of the electron 
fluid near the boundary of the gated region. The spec- 
trum of the plasma oscillations should be broad, similar 
to what is observed in experiments (Fig. 1). The width 
of the spectrum is expected to be limited by the value 
Umax ~ s/d, where d is the gate-to-channel separation 
(see above). 

The present theory can be also applied to the ungated 
electron fluid (analogous to the "deep water" in conven- 
tional hydrodynamics). It was shown [^] that a one- 
dimensional instability similar to the one described in 
Ref. [l| should exist in the ungated case too, under ap- 
propriate boundary conditions. It can be easily shown, 
that the boundary instability considered here will also 
occur in the ungated region, if the drift velocity is di- 
rected inside this region, similar to the results given by 
Eqs. (17,20) for the gated electron fluid. Thus, at the 
boundary between the gated and ungated regions, the 
turbulence should always appear on the downstream side. 

It appears that the above concept acounts for the most 
important experimental observations [H, 0, [1[ : the sharp 
threshold for terahertz emission and the broad emission 
spectrum, which does not depend on the gate length, and 
only weakly depends on the gate potential. A possible 
check of the proposed explanation would be to isolate 
the emission coming from one gate edge and to verify 
that the emission intensity (and possibly its spectrum) 
depends on the direction of the the drift velocity. 

On the theoretical side, the very difficult issue of the 
true conditions at the boundary between the ungated and 
gated regions should be elucidated. (From the hydrody- 
namical point of view this is the problem of what happens 
for a flow across the boundary between deep and shallow 
water). Also, the role of the viscosity of the electron fluid, 
which may supress the instability for large wave vectors 
q, remains to be understood. 
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